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具有 质 体 和 抑制 剂 的 非 均匀 恒 化 器 模型 的 分 歧 * 


WRI, x UOS, REH? 
(1 西安 石油 大 学 理学 院 ， 西 安 710065; 2- 陕西 师范 大 学 数学 与 信息 科学 学 院 ， 西 安 710062) 


fü 要 : 本 文 讨 论 了 一 类 具有 外 加 抑制 剂 的 质 体 负载 与 质 体 自由 物种 竞争 的 非 均 匀 恒 化 器 模型 共存 解 的 存 
在 性 和 稳定 性 。 通 过 比较 原理 、 分 歧 理 论 和 线性 稳定 性 理论 ， 分 析 了 由 半 平 凡 解 发 出 的 分 歧 解 的 
全 局 结构 和 局 部 稳定 性 ， 解 释 了 物种 共存 的 现象 。 结 果 表 明 抑 制剂 有 利于 遗传 选择 的 物种 避免 被 


质 体 自由 的 物种 所 消亡 。 
关键 词 : 恒 化 器 ; 质 体 负 和 载 ， 质 体 自由 ; 共存 解 ， 全 局 分 歧 
分 类 号 : AMS(2000) 35B40 中 图 分 类 号 : O175.26 文献 标识 码 : A 
1 引言 


近来 ， 具 有 质 体 的 恒 化 器 竞争 模型 受到 了 广泛 的 关注 [四 。 本 文 研究 如 下 具有 质 体 和 外 加 抑 
制剂 的 非 均 匀 恒 化 器 竞争 模型 
St = dsa4 — laufi(s) — lhvfa(s)e-"?, r€(01) t0, 
ut = dus + a(1— q)ufi(s), z€(0,.1), t0, 
vt = duza + bvfa(s)e "? + aqufi(s), z€(01) #1>0, 
pi = dpzz — cuh(p), rE€(0,1), Ł>0, 


(1) 


具有 边界 条 件 


sz(0,t) = 一 s0， Sz(1,t) + ys(1,.) 20, t>0, 
p. (0,t) = —p, pz(1,t) - yp(1,t) 20, £2 0, (2) 
us (0, t) = us (1,t) + yu(1, t) = 0, v; (0, t) = v. (1,t) + vv(1,t) —0, t»0 


和 非 负 的 初始 条 件 。 这 里 (m, 0), u(x, t), v(x, t), p(x,t) 分 别 为 :时刻 营 养 物 、 两 竞争 物种 和 撮 
制剂 的 浓度 。s0, p9 > 0 均 为 常数 ， 分 别 为 营养 物 和 抑制 剂 的 输入 浓度 。d 为 恒 化 器 的 扩散 
E, 7 为 产 出 率 。a,b 为 两 竞争 物种 的 最 大 生长 率 ， 反 应 函数 fils) = s/(ai + s)(i = 1,2), K 
中 ai 为 半 饮 和 常数 ，c 为 物种 对 换 制 剂 的 吸收 率 ，h(p) = p/(hi 十 p)， 其 中 hi 为 半 饱 和 常 
数 。 记 质 体 游离 的 概率 为 gs， 则 0 < gq < 1. e-^? 体现 了 抑制 剂 z 对 物种 v 生 长 率 的 抑制 程度 ， 
Hh u> 0 为 常数 ， 反 应 了 抑制 剂 对 物种 v 的 影响 。 

基于 文献 3,4] 的 研究 结果 ， 许 世 壁 等 丫 提出 了 一 类 具有 质 体 和 外 加 抑制 剂 的 均匀 搅拌 的 恒 
化 器 模型 ， 得 到 了 此 模型 的 一 些 渐 近 行 为 。 文 献上 改进 了 他 们 的 结果 ， 得 到 了 模型 存在 唯一 


收 稿 日 期 : 2008-11-10. 作者 简介 : 谢 文 丘 (1978 年 10 月 生 )， 女 ， 硕 士 ， 讲 师 . 研究 方向 ， 偏 微分 方程 数值 解 . 
* 基 金 项 目 : 国家 自然 科学 基金 (10971124); 陕西 省 自然 科学 基础 研究 计划 (2009JQ1007). 


第 5 期 谢 文 吴 等 ， 具 有 质 体 和 抑制 剂 的 非 均匀 恒 化 器 模型 的 分 歧 839 








正平 衡 解 的 充 要 条 件 以 及 模型 存在 周期 解 的 充分 条 件 。 文 献 [6] 去 掉 均匀 搅拌 的 假设 ， 考 虑 了 
一 类 具有 质 体 和 内 部 抑制 剂 的 非 均匀 恒 化 器 模型 ， 给 出 了 当 抑制 剂 的 影响 足够 大 时 ，PDE 系 
统 的 一 些 渐 近 性 。 本 文 考虑 外 加 抑制 剂 的 情形 ， 即 研究 如 下 简化 的 椭圆 系统 


du” 4 a(1— q)ufi(z — u — v) — 0, u'(0) = w'(1) + yu(1) = 0, 
dv" + bvfa(z — u — v)e-"? c aqufi(z —u— v) 20, v'(0) — v'(1) + yv(1) — 0, (3) 
dp. — cuh(p) = 0, p(0)— —1, p'(1)-* vr(1) = 0, 


共存 解 的 分 支 结 构 和 稳定 性 ， 其 中 z(z) = 1 — z。 这 里 (3) 是 在 (1)- (2) 的 平衡 态 系统 中 引入 
无 量 纲 参数 和 变量 ， 并 使 用 恒 化 器 模型 所 固有 的 守恒 率 而 得 到 的 ， 详 细 的 推导 见 文献 [5-7]。 

由 于 我 们 仅 关 心 (3) 的 共存 态 ( 即 (3) 的 正解 )， 因 此 可 将 反应 函数 延 拓 为 C1 连续 的 函数 回 。 
下 面 引 入 一 些 记 号 和 一 些 已 知 的 结论 。 设 Xi m 分 别 为 如 下 特征 值 问 题 的 主 特征 值 





de, + Afilz)p1 ==0，ze(0,1)， pi (0) = v1) + yp1(1) = 0, 
di, - mfa(z)e "^W. —0, z€(01) (0) =v (1) - vi (1) — 0, 


相应 的 正 单 位 特征 函数 满足 maxi — i; meyi = 1。 考 虑 如 下 退化 的 非 线性 边 值 问题 


du -a(1l—q)ufi(z—u)-0, z€(01) w(0)—w(1)4vu(1) — 0, (4) 

dv -bvofa(z—v)e""* 20, z€(01)  v'(0) =v (1) + yv(1) — 0. (5) 
显然 ， 若 a < X/(1—q), Wu = 0 为 (4 的 唯一 非 负 解 ; 若 w A/(1— a), W (4) 存在 唯一 正 
解 ， 记 为 Wa)， 且 0 < 3(a) < z (参见 文献 [7])。 而 且 


£x) id 
在 [0,1] 上 一 致 成 立 。 类 似 地 ， 若 b xo. Bv = 0 为 (5) 的 唯一 非 负 解 ; ib > m. BU (5) 存在 
唯一 正解 6， 且 0 < 0 < z。 

本 文 主要 研究 系统 (3) 共存 态 的 全 局 结构 和 稳定 性 。 由 于 0 < g < 1， 即 质 体 可 能 游离 ， 使 
有 具有 质 体 (plasmid-bearing) 的 物种 转化 为 不 具有 质 体 (plasmid-free) 的 物种 ， 这 就 破坏 了 标准 
的 恒 化 器 模型 所 具有 的 严格 竞争 性 ， 因 此 单调 系统 理论 在 这 里 失效 。 但 借助 于 一 些 重要 的 特征 
值 估计 和 比较 原理 ， 仍 然 可 以 证 明 上 > n b < m 时 (3) 均 存在 共存 解 。 


2 ”共存 解 的 存在 性 与 稳定 性 


本 节 主 要 讨论 系统 (3) 共存 态 的 全 局 结构 和 稳定 性 。 为 此 ， 首 先 讨论 (3) 正解 的 先 验 估计 。 
引 理 2.1 如 下 边 值 问题 存在 唯一 正解 pa。， 且 0 < pa < zo 


dp —có(a)h(p —0, z€(01)  p»(0--1 p(1)+7yp(1)=0. (6) 


证 明 设 p(z) 为 (6) 的 非 负 解 ， 则 可 以 证 明 0 < p(z) < z(z)。 实 际 上 ， 由 最 大 值 原理 
有 p(x) > 0(z € [0,1)。 令 g= z(7x) 一 p， 则 4 满足 


-dd = cd(a)h(p) > 0, q'(0)-— q'(1) + ya(1) = 0. 
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根据 最 大 值 原理 可 知 g(z) > 0(z € [0,1)， 即 0 < plz) < z(z)。 另 一 方面 ，z(z),0 分 别 为 (6) 的 
上 下 解 。 因 此 (6) 存在 正解 z(z)， 且 0 < p < z。 下 面 证 明 (6) 在 (0,z(z)) 内 只 有 一 个 正解 。 
令 0< pl(z)，pa(z) < z 为 (6) 的 两 个 正解 ， 则 g = pi — p 满足 


" chi(a) 
7 Ts ph tp)! 


再 由 最 大 值 原 理 有 q = 0， 即 唯一 性 成 立 。 
下 面 给 出 (3) 正解 的 先 验 估计 ， 其 证 明 参 见 文献 [7] 的 引 理 4.1 和 引 理 4.2。 
引 理 2.2 设 (w up) 为 (3) 的 非 负 解 且 满足 #0, v 0, W 
(i) a>A/(l-g); (Gii) 0<u< Va), 0cv«z,p,«p-«zi: (üi) 0<u+v< zo 
引 理 2.3/5! ibm. UU 


0, z€(01) g'(0) =q (1) + ya(1) = 0. 


d? 
Le = d 十 b (falz d 8) = 8 falz E 0)).e"* 


的 所 有 特征 值 都 是 负 的 。 
下 面 采 用 分 歧 理 论 研究 (3) 正解 的 结构 。 设 


C$ (I0, 1]) = {u € C7{([0,1]) : v'(0) = v'(1) + Yu(1) = 0) 


有 具有 通常 的 C- A |- l, X = C801) x C5(0.1) x CB([0, 如 )。 首 先 考虑 b > o 的 情形 ， 
于 是 系统 (3) 存 在 平凡 的 非 负 解 (0,0, z) 和 半 平 凡 的 非 负 解 (0,9,z)。 取 a 为 分 歧 参 数 ， 从 半 
平凡 的 非 负 解 分 支 {(a,0,0,z) : a € Rt+} 出 发 构造 (3) 的 正解 。 令 w= wu, x=0-v, ¿= 
z—p KJ -d fECL(0,1]) PREET. WF = (F, Fa Fa) EP | 


Fi(w; x,&) 7 a(1— ajo(fi(z — 9 — w + x) — fi( — 0)), 
Folw, x, £) = b(x — 0)fa(z — 0 — w + x)e "979 — b(x — 0) falz — Oe 
+ b6f5(z — Oe "(x — w) + aqw(fi(z — 8 -w + x) — fi(z — 0)) 
+ pbb fo(z — 0)e "^£, 
Fa(w, x, £) = ew(h(z — £) — h(2)), 
N FÆ, F(0,0,0)—0, H Fréchet 导数 Dc, ,F(0,0,0) 20. SIET, GA 
a(1 — q)K (wfi(z — 0)) + KFi (w, x, £), 
bK (x(fa(z — 0) — 8f2(z — 6))e"*) + K[(bOfs(z — 0)e-"* 
T(a,w,x,£) = 
— agfi(z — 0))w] — ubK (8f2(z — 0)e-"*£) + K Fo(w, x, €), 
cK (h(z)w) + KFs(w,x,€) 
G(a,w, x £) = (w, x: £) m T(a,w, x: £), 


则 工 (a,w, xX,é&) 为 XX 上 的 连续 可 微 紧 算 子 ，G 是 C1 连续 的 ，G(a,0,0,0) —0. ifi H.G(aw,x,£) 
二 0 满足 0 <w < Va), 0<0 一 x <z, 0<z--&<z 的 解 恰好 为 (3) 的 正解 。 
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考虑 在 点 (a.u, x, E) = (E 1,0,0,0) 的 分 歧 ， 其 中 À 为 如 下 边 值 问题 的 主 特征 值 
d$; F Afi(z = 0)P1 = 0, TE (0, 1), (0) = ĝi (1) T YP1(1) 一 0， 
相应 的 主 特征 函数 2 > 0 (zx € [0,1)， H maxĝ, = =1.  L(a,0,0,0) = Dwxe)G(a,0,0,0) 为 
G XT (cw, x, €) 4€ (a,0,0,0) ff] Fréchet 导数 ， 则 易 验证 核 空间 


N(L( 0, 0, 0)) = span((à1, xi, &)), 














值 域 & T 
A(Z( -2 31000))= (e «nex: f unas - 0), 
mH 
L( 0, 0, 0) ($1,231, &) = Dc, je es ; 0, 0, 0) (£1, xi, &) 
- es 0)), qK ($1f1(z — 0)),0) 
- (- X0 - 96 Ao 0) ¢ n(1( 72, ,0,0) )， 
其 中 


&= K (ch(z)ĝ1) >0, 





vsi Du file - 9i) + Ez (9(— Fle- 9)  nfaz — 06)e7"" ]. 
由 文献 [8] 的 定理 13.5 TAI RA (335,0, 0, z) 的 (3) 的 正解 分 支 为 光滑 曲线 


Fi = ((a(s), s($1 + 9(s)), 0 — soa + v(s)), z — s(& + Q(s))) : s € (0,8)}, 











«(0)- 3, 4()-9()-Q0)-0, (49) W), Q(9) € {Cr x £)]- 
进一步 ， 类 似 于 文献 6, 可 以 证 明 a(s) 关于 s 在 s = 0 点 的 导数 为 


Ado f(z - 0): - d (n) 
sed 1-9); filz- 
下 面 证 明 Ti 可 延 拓 为 全 局 分 支 。 设 ai( 入 ) 为 如 下 问题 的 特征 值 
Adé taNfi(z~060=0, ze(0,1), d(0)—dQ(1) -y9(1) =0, 
则 ailà) 在 (1, +oo) 内 关于 入 严格 递增 , EURO < ali(A) < az(A) < as(A) < … 一 œ, a(1) = 
Al。 根据 基本 的 全 局 分 歧 理 论 67] 可 以 证 明 当 a < A NS i(T(a,),0) = 1; 当 


da 
ds 





à aj(1) 
jog 0s TEN 


8 


Bf, d(T(a,.),0) = -1。 由 文献 们 的 全 局 分 歧 定理 2.1 可 知 ， 在 Rx X 内 存在 发 自 ( 疡 - ,0,9,z) 的 
连续 分 支 满 足以 下 三 种 情形 之 一 ， 

(i) 工 连接 了 点 (2 0,b,z)， 其 中 五 天 ju; 

(i) TÆR+ x X Ah (2,0,0, 2) 延伸 到 无 穷 

(ui) D &SSJÉUI(a,u,0 — v, z — p) fll (a, —u,0 +v,z +p) WA, HF (u,v, p) 去 (0,0,0). 


382.4 T- ((225,0,6,2)) Hi (235,0,0, 2) 出 发 在 正 锥 忆 内 延伸 到 无 穷 ， 其 中 


P = {(a,u,v,p)E Rt xX:u>0,v>0,p>0,ze [0,1]}. 


证 明 假设 
A 
r-[(12566:))g» 
则 存在 . 
VE, m 
(à,à, dD (T — (G259*.2))) NaP 
MIES {Ön Ün Ün Pn)} C TO PP 使 得 当 n — oo. (nün Ün Pn) — (ü, p) Bc HU 
大 值 原 理 可 知 方 0. AE (á, ü, 0,5) € 9P 表 明 存 在 zo € [0,1] Ei > 0, ü(zo) = 0, 


ü 
或 者 存在 zo € [01] 使 得 5 > 0, (zo) = 0。 再 根据 最 大 值 原理 有 立 三 0 或 者 5 = 0。 显 
然 ,， 0—0, 让 > 0 不 成 立 。 

假设 (ü, V) 三 (0, 0); Ju 35 n = co 时 ， (bn, ün, Ün, Pn) E (à, 0, 0, z)e 而 且 Ün = iin / | s ||, 
满足 


dU, -à,(1—g)Ü,fi(z— à, —9,)—0, z€(01, 0UL(0) = Č}(1) + 0,(1) = 0. 
H L? 估计 和 Sobolev KAEST AI, FE u € CL(0,1]) Wig u 0, uz, H 
du -à(l—q)ufi(z)—0, z€(0,1) w(0) =u (1) + yu(1) = 0. 
根据 最 大 值 原理 有 > 0 (x e [0, 1])， 这 表明 & = 入 1/(1 一 gq)。 于 是 
à, — X/(0 —q)(n— oo), lim Hn) = 0, 


其 中 (G64) 为 a = à, RT (4) 的 唯一 正解 。 
由 引 理 2.2 可 知 


di;, 十 bó, falz —- Ün P Ün )e P» T ü, qi, fi(z t Ün Ee Ün) > di,, T bb, f2 (z -— Ü(à,.) Er. Ün jez. 
因此 ， 如 果 > M(fa(z—9(a,)e-"*), MG, > wn > 0， 其 中 性 i 为 
dà, T bo, f2(z = 2(Gn) e is )B. A7 = 0, tE (0, 1), wn(0) = wn.(1) 十 Yon(1) =0 


的 唯一 正解 。 实际 上 和 (fo(z — 9(à5))e^"?) — m(n 一 oo)。 X bg 因此 存在 入 > 0， 使 得 
Zn2NBbA(f(z-39(u4)e ). FEÆïn > Nitin 0,0. X 9, 一 0(n — oo), 
Eb o, 一 0(n oo). S Wn —o./|ios]. JU Wn WE 


dW, --bW,fo(z — din) — õn) "* —0, z€(01) WZ(0-W.)-44W,(1)- 0. 
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1848 L^ hiit A Sobolev RA EM, Ew > 0, w z 0, EE CLIO, 1]) HW, 一 win 一 
œ) H 
dw” +bwfa(z)e™™* 20, ze (0,1), w'(0) =w (1) +yw(1) = 0. 


而 且 ， 由 最 大 值 原理 可 知 w > 0(z e [0,1])， 这 表明 b= 二 1， 与 b> mJ. 
假设 让 三 0, 证 > 0， 则 = 90, p= zo Bl (än, ün, 95, Du) > (6,0,0,2) (n — oo)» $Ü, = 
ü,/|ü,], WU, 满足 


dÜ, rü,(1—g)Ü,fi(z— à, —9,)—0, z€(01), U/(0 = ČL(1) + 40,(1) = 0 


同 理 可 证 & = Èy 了 矛盾。 

TJ&r-((22,0,0,2)) C P, TUE P KERIEN. MEHUR (u,v, p) ET-{(ao,0,0,z)}， 
HBI2.2, LP 估计 和 Sobolev 嵌 入 定理 ， 则 存在 M > 01558 lull, loll, lip] < M。 因 此 ， 在 正 
HE PAT UW a REN. 

综 上 ， 我 们 有 如 下 全 局 分 歧 的 结 

定理 2.1 b > m， T A 5,0,0,2) 的 (3) 正 解 的 连续 分 支 P， 它 在 正 锥 尸 内 
Wa 延伸 到 无 穷 ， 且 Ti Cr. 

下 面 研 究 分 歧 点 (这 pho 2) 附近 共存 解 的 稳定 性 。 令 Xi = [C?^((0,1)? n X, Y = 
[C*([0, 1)]?, i: X1 一 工 为 包含 映射 。 类 似 于 文献 四] 的 引 理 5.2 易 知 ， 029 L( 3-0, 0,0) 85 i- 
RRE EL L(,315,0,0,0) 的 所 有 其 它 特 征 值 都 位 于 左 半 复 平面 。 

V L(a,0,0, z) Wl L(a(s), u(s), v(s), p(s)) 分 别 为 (3) 在 点 (ao,0,0,z) II (a(s), u(s), v(s), p(s)) 
处 的 线性 化 算 子 。 根 据 文献 [8] 的 引 理 13.7 和 引 理 13.8 可 知 ， 存 在 分 别 定义 在 A ; 和 0 邻 域 上 
的 C1 连续 函数 a 一 (n(a), 9(a)) e Rx Xi 和 ss 一 (r(s), (s) e Rx Xy, 1% 


(7( e( ^ )) = (0, à, X1, &) = (7(0), w(0)), 





ic 
H a 
L(a,0,0, z)*(a) = v(a)9(a), a 一 zl « 1, (8) 
L(a(s), u(s), v(s), 2D(s)) 亚 (s) = r(s)W(s) |s| «1, 
其 中 


(a) = ($:(a), 92(a), Ba(a)), 亚 (s) = (Wi(s), V2(s), Vs(s)). 

idm j £90 X4r(s) # 0(|s| « 1)， 则 当 |s| « 1 时 7(s) 和 一 -sà(sji( È) 同 号 ， 其 

Hre JA rla) X Fatta = 地 处 的 导数 ， &(s) 为 a(s) 关 于 s 的 导数 。 进 一 步 ， 类 似 于 文 
献 [5] 的 引 理 3.7 可 以 证 明 方 (aa 5) > 0。 再 结合 等 式 (7)， 有 如 下 稳定 性 定理 。 

定理 2.2 令 
I= 'ga 1(z — 9)(91 — x1)d 
- f eif Q1-— Xi 。 

$120, Wjr(s) <0(0 <s < 1)， 因 而 定理 2.1 所 给 出 的 分 歧 解 稳定 。 

注 2.1 HURDA, HIULUEB 0I Ly bol f2(z — 0)1 H+ ufa(z — 0)61)e-"?] JF 5E, 
HFE ô > 0 使 得 





入 
一 1 1q " 
万 (Tf -0 ) « 6. 
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又 因为 & = K(ch(z)$1) > 0， 因 此 只 要 1 充分 大 ， 就 有 Xi < 0， 从 而 有 7 > 0， 且 由 定 
理 2.1 给 出 的 分 歧 解 在 0 < s < 1 时 稳定 。 从 生态 学 上 看 ， 参 数 体现 了 抑制 剂 的 影响 。 因 此 这 
一 结果 正好 表明 抑制 剂 有 利于 遗传 选择 的 物种 避免 被 不 具有 质 体 的 物种 所 消亡 。 

对 于 b < 六 的 情形 ， 同 理 有 如 下 结论 。 

定理 2.3 设 b < m JU (3) 存在 发 自 ( 产 -,0,0,z) 的 正 全 局 分 支 F， 它 在 已 内 沿 a 延 伸 到 
无 穷 。 而 且 ， 在 分 歧 点 ( 疡 5,0,0,z) 附近 ， 正 解 稳定 ， 且 


A 


(a: (a,u, v, p) E} = ert o). 


注 2.2 定理 2.3 表 明 如 果 5 < 六 固定 ， 则 a € (u-, o0) 时 (3) 存在 共存 解 。 


1 一 和? 
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Bifurcation on the Plasmid-bearing and Plasmid-free Model in the 
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Abstract: The existence and stability of coexistence solutions to a competition model between 
plasmid-bearing and plasmid-free organism in the unstirred chemostat with an external inhibitor are 
discussed. By the comparison principle, bifurcation theory and linear stability method, the global 
structure of the bifurcation solution from semitrivial solution and their local stability are analyzed, 
which explain the phenomenon of coexistence. The results indicate that the inhibitor can help the 
genetically altered (plasmid-bearing) organism to avoid 'capture' of the process by the plasmid-free 
organism. 
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